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' Ground states of four-dimensional (d=4) EA Ising spin 

glasses are calculated for sizes up to 7x7x7x7 using a 
combination of a genetic algorithm and cluster-exact approx- 
imation. The ground-state energy of the infinite system is 
extrapolated as eo° ~ —2.095(1). The ground-state stiffness 
(or domain wall) energy A is calculated. A |A| ~ L®^ be- 
havior with Qs = 0.65(4) is found which confirms that the 
d=4 model has an equilibrium spin-glass-paramagnet transi- 
tion for non-zero Tc. 

Keywords (PACS-codes): Spin glasses and other 
random models (75.10.Nr), Numerical simulation studies 
(75.40.Mg), General mathematical systems (02.10.Jf). 



I. INTRODUCTION 

Optimization methods have found widespread applica- 
tion in computational physics. Among these the inves- 
tigation of the low-temperature behavior of spin glasses 
attracted most of the attention within the statistical 
physics community. The reason is that despite its sim- 
ple definition (see below) its behavior is far from being 
understood. From the computational point of view the 
calculation of spin-glass ground states is very demanding, 
because it belongs to the class of the NP-hard problems 
This means that only algorithms are available, for 
which the running time on a computer increases expo- 
nentially with the system size. In this work a method re- 
cently proposed, the cluster-exact approximation (CEA) 
1^ is applied to four-dimensional Ising spin glasses. 

The model under investigation here consists of N spins 
CTi = ±1, described by the Hamiltonian 



H 
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where (...) denotes a sum over pair of nearest neigh- 
bors. In this report simple 4d lattices are considered, i.e. 
TV = L"*. The nearest neighbor interactions (bonds) take 
independently = ±1 with equal probability. Periodic 
boundary conditions are applied to the systems. No kind 
of external magnetic field is present here. 

Four-dimensional Ising spin glasses have been investi- 
gated rather rarely. Most of the results were obtained via 
Monte-Carlo (MC) simulations at finite temperature, see 
e.g. P-p^. Here the T = behavior is investigated, i.e. 
ground states are calculated. This has the advantage. 



that one does not encounter ergodicity problems or crit- 
ical slowing down like in algorithms which base on MC 
methods. Only one attempt to address the 4d spin- 
glass ground-state problem is known to the author. But, 
as we will see later, the former results suffer from the 
problem, that not the true global minima of the energy 
were obtained. Furthermore, no analytic predictions of 
the ground-state energy have been noted by the author. 

The question whether finite-dimensional Ising spin 
glasses show an ordered phase below a non-zero transi- 
tion temperature Tc is of crucial interest. By MC simula- 
tions around the (expected) transition temperature this 
question is hard to solve. Another way to address this 
question is to calculate the stijfness or domain wall en- 
ergy A — i?" — EP which is the difference between the 
ground-state energies E°',EP for antiperiodic and peri- 
odic boundary conditions in one direction Here 
the antiperiodic boundary conditions for calculating 
are realized by inverting one plane of bonds. For the 
other directions periodic boundary conditions are applied 
always. This treatment introduces a domain wall into the 
system. If a model exhibits an ordered low-temperature 
phase, the domain wall increases with growing system 
size, which becomes visible through the behavior of A: 
the disorder-averaged stiffness energy shows a finite size 
dependence 
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A positive value of the stiffness exponent Qs indicates 
the existence of an ordered phase for non-zero tempera- 
ture. For example a simple d = 2 Ising ferromagnet has 
O5 = 1. For spin glasses, the stiffness exponent plays 
additionally an important role within the droplet-scaling 
theory [pf-p^ , where it describes the finite-size behavior 
of the basic excitations (the droplets). 

Using this kind of analysis is was proven that the 2d 
spin glass exhibits no ordering for T > |1^]. For the 
three-dimensional problem in a recent calcu ation ]20| by 
applying genetic CEA a value of Og ~ 0.19(2) was found, 
which shows, that indeed the d = 3 model has a spin-glass 
phase for nonzero temperature. For d = A the existence 
of a finite Tc ~ 2.1 was proven rather early even by MC 
simulations but the value for the stiffness-exponent 
Qs is of interest on its own. In ||l^ recently a value of 
Qs = 0.82(6) was found by performing a MC simulation 
near Tc- In the work presented here the value is obtained 
via ground-state calculations. 
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The paper is organized as follows: In the next section 
the algorithm applied here is briefly presented. The main 
section contains the results for the ground-state energy 
and the stiffness exponent. Finally a summary is given. 
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III. RESULTS 



II. ALGORITHM 

The technique for the calculation bases on a special 
genetic algorithm ||2l|,^ and on cluster-exact approxi- 
mation 1^ which is an optimization method designed es- 
pecially for spin glasses. Now a brief description of the 
method is given. 

Genetic algorithms are biologically motivated. An op- 
timal solution is found by treating many instances of the 
problem in parallel, keeping only better instances and 
replacing bad ones by new ones (survival of the fittest). 
The genetic algorithm starts with an initial population of 
Mi randomly initialized spin configurations (= individu- 
als), which are linearly arranged in a ring. Then v x Mi 
times two neighbors from the population arc taken (called 
parents) and two offspring are created using the so called 
triadic crossover |p3|| . Then a mutation with a rate of pm 
is applied to each offspring, i.e. a fraction pm of the spins 
is reversed. 

Next for both offspring the energy is reduced by ap- 
plying CEA. The algorithm bases on the concept of frus- 
tration |2^. The method constructs iteratively and ran- 
domly a non-frustrated cluster of spins, whereas spins 
with many unsatisfied bonds are more likely to be added 
to the cluster. The non-cluster spins act like local mag- 
netic fields on the cluster spins. For the spins of the clus- 
ter an energetic minimum state can be calculated in poly- 
nomial time by using graph-theoretical methods p5|-^7[ : 
an equivalent network is constructed [p8[|, the maximum 
fiow is calculated ^9 30 1 and the spins of the cluster are 
set to the orientations leading to a minimum in energy. 
This minimization step is performed n^i^ times for each 
offspring. 

Afterwards each offspring is compared with one of its 
parents. The pairs are chosen in the way that the sum 
of the phenotypic differences between them is minimal. 
The phenotypic difference is defined here as the number 
of spin where the two configurations differ. Each parent 
is replaced if its energy is not lower (i.e. better) than the 
corresponding offspring. 

After this creation of offspring is performed h' x Mi 
times the population is halved: From each pair of neigh- 
bors the configuration which has the higher energy is 
eliminated. If not more than 4 individuals remain the 
process is stopped and the best individual is taken as 
result of the calculation. 

The whole algorithm is performed un times and all 
configurations which exhibit the lowest energy are stored, 
resulting in rig statistical independent ground state con- 
figurations. The method was already applied for the in- 
vestigation of the ground-state landscape of 3d Ising spin 



In this section at first the values for the simulation pa- 
rameters, which are defined above, are presented. Then 
the finite-size behavior of the ground-state energy is in- 
vestigated. Finally results for the stiffness energy are 
discussed. 

The simulation parameters were determined in the fol- 
lowing way: For the system sizes L — 2,4,6,7 several 
different combinations of the parameters Mi,v,nmin,Pm 
were tested. For the final parameter sets it is not possible 
to obtain lower energies even by using parameters where 
the calculation consumes four times the computational 
effort. For L — 2>,b the parameter sets for i + 1 were 
used. Using parameter sets chosen this way genetic CEA 
calculates true ground states, as shown in |2y]. It should 
be pointed out that it is relatively easy to obtain states, 
which exhibit an energy slightly above the true ground 
state energy. The hard task is to obtain really the global 
minimum of the energy. 

Here Pm = 0.1 and n/j = 5 were used for all system 
sizes. Table Q summarizes the parameters. Also the typ- 
ical computer time r per ground state computation on a 
80 MHz PPC601 is given. 

Ground states were calculated for system sizes up to 
7 X 7 X 7 X 7 for Nl independent realizations (see table |) 
of the random variables. For each realization the ground 
states with periodic and antiperiodic boundary condition 
in one direction were calculated. The remaining three di- 
rections are always subjected to periodic boundary con- 
ditions. One can extract from the table that for small 
system sizes L < A ground states are rather easily to ob- 
tain, while the L — 7 systems alone required 6560 CPU- 
days. Using these parameters on average rig > 2.7 ground 
states were obtained for every system size L using nn = 5 
runs per realization. 

The average ground-state energy cq per spin is shown 
in Fig.|l| as a function of the system size L. Using a 
fit to eo{L) ~ + a * L^^ the value for the infinite 
system is extrapolated, resulting in = —2.095(1) 
(a = 7.1(7), 5 = -4.2(1)). This value is compatible with 
the lower bound of ep — \/2d\n2 sa 2.35 given by the 
random energy model [Q. The value calculated here is 
substantially smaller than the result eg° = —2.054(3), 
which was obtained in pH using a pure genetic algo- 
rithm. This shows that in ||ll| not the true global min- 
ima were found, which can be concluded also from the 
fact, that there eo{L) increases with growing system size. 
Because the periodic boundary conditions impose addi- 
tional constraints on the systems, the opposite behavior 
is expected, as found for the results presented here. For 
further comparison additionally some calculations were 
performed by the author by simply rapidly quenching 
from random chosen spin configurations. By executing 



2 



an analogous fit, a value of eg° = —2.04(2) is obtained. 
This shows, that the result from |0 seems to be only 
slightly better than the data obtained by applying a very 
simple minimization method. 

The distribution of the stiffness energy, which is ob- 
tained from performing ground-state calculations for sys- 
tems with either periodic or antiperiodic boundary con- 
ditions in one direction, are shown in Fig. |^ for L — 5 
and L = 7. With increasing system size the distribution 
broadens. This means that larger domain walls become 
more and more likely. To study this effect more quanti- 
tatively, in Fig. ^ the disorder-averaged absolute value 
(|A|) of the stiffness energy is plotted as a function of the 
system size L. Also shown is a fit (|A(L)|) ~ i®^ which 
results in 8s — 0.64(5). Here, the system sizes L — 2,3 
were left out of the analysis, since they are below the 
scaling regime. Because of the large sample sizes the er- 
ror bars are small enough, so we can be pretty sure that 
05 > 0. It confirms earlier results from MC simulations 
|,| that the 4d EA spin glass exhibits a non-zero transi- 
tion temperature Tc- The value Qs = 0.64(5) is compara- 
ble to a recent result from MC simulations 85 ~ 0.82(6) 
pof , given the facts that the system sizes are rather small 
and the other result was obtained at finite temperature 
near the transition point Tc ~ 2.1. Additionally, the pre- 
diction from droplet-scaling theory 8s < (d — l)/2 = 1.5 
is fulfilled. 

It should be pointed out, that the method described 
above does not guaranty to find exact ground states, al- 
though the method for choosing the parameters makes 
it very likely. If states with a slightly higher energy are 
obtained, the result for ejf is not affected very much. For 
the stiffness energy, it was shown in |2^ that the result is 
very reliable as well, as long as the energies of the states 
are not too far away from the true ground-state energies. 



IV. CONCLUSION 

Results have been presented from calculations of a 
large number of ground states of 4d Ising spin glasses. 
They were obtained using a combination of cluster-exact 
approximation and a genetic algorithm. Using a huge 
computational effort it was ensured that true ground 
states have been obtained with a high probability. 

The finite size behavior of the ground-state energy and 
the stiffness energy have been investigated. By perform- 
ing a L 00 extrapolation, the ground-state energy 
per spin for the infinite system is estimated to be = 
—2.095(1). The absolute value of the stiffness energy in- 
creases with system size and shows a (|A(L)|) ~ L®^ be- 
havior with 8s — 0.64(5). For systems with a Gaussian 
distribution of the bonds qualitatively similar results are 
expected, since the ordering behavior depends only on 
the sign of the interactions and not on their magnitudes. 

A more detailed study of the ground-state landscape 
of 4d systems, similar to pl|, requires more than no ~ 



3 ground states per realization to be calculated. Since 
this requires a substantial higher computational effort, it 
remains to be done for the future. 
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FIG. 2. Distribution of the stiffness energy A = - 
for system sizes 5x5x5x5 and 7x7x7x7. E" and E'^ are 
the total ground-state energies for periodic and antiperiodic 
boundary conditions in one direction, while for the other three 
directions always periodic boundary conditions are imposed. 
Lines are guide to the eyes only. 
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FIG. 1. Average ground-state energy eo per spin as 
a function of system size L. The line shows a fit to 
eo{L) = +a*L~'' resulting in eg° = —2.095(1) as estimate 
for the ground-state energy of the infinite system. 
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fit: <|A|>=aL''s 
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©3=0.65(4) 
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FIG. 3. Average Stiffness energy (|A|) as function of sys- 
tem size L on log-log scale. The line represents the function 
|A(L)| = aL®s with Os = 0.65(4). The increase of (|A|) 
with system size indicates, that for 4d Ising spin glasses an 
ordered phase exists below a non-zero temperature T^. 
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TABLE I. Simulation parameters: L = system size, M, = 
initial size of population, v = average number of offspring per 
configuration, rimin = number of CEA minimization steps per 
offspring, t = typical computer time per ground state on a 
80MHz PPC601, Nl = number of realizations of the random 
variables. 
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